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Abstract 

We consider approximate cloaking from a regularization viewpoint introduced 
in [13] for EIT and further investigated in T2l [17] for the Helmholtz equation. 
The cloaking schemes in [12] and [17 are shown to be (optimally) within | In p\~^ 
in 2D and p in 3D of perfect cloaking, where p denotes the regularization param- 
eter. In this paper, we show that by employing a sound-hard layer right outside 
the cloaked region, one could (optimally) achieve in M.^ , N > 2, which signif- 
icantly enhances the near-cloak. We then develop a cloaking scheme by making 
use of a lossy layer with well-chosen parameters. The lossy-layer cloaking scheme 
is shown to possess the same cloaking performance as the one with a sound-hard 
layer. Moreover, it is shown that the lossy layer could be taken as a finite re- 
alization of the sound-hard layer. Numerical experiments are also presented to 
assess the cloaking performances of all the cloaking schemes for comparisons. 

1 Introduction 

A region is said to be cloaked if its contents together with the cloak are invisible to 
certain measurements. From a practical viewpoint, these measurements are made in 
the exterior of the cloak. Blueprints for making objects invisible to electromagnetic 
waves were proposed by Pendry et al. [23] and Leonhardt |16j in 2006. In the case 
of electrostatics, the same idea was discussed by Greenleaf et al. PjDj in 2003. The 
key ingredient is that optical parameters have transformation properties and could 
be pushed- forward to form new material parameters. The obtained materials/media 
are called transformation media. We refer to [H [3 El [231 [261 EE] for state-of-the- 
art surveys on the rapidly growing literature and many striking applications of the 
so-called 'transformation optics'. 

In this work, we shall be mainly concerned with the acoustic cloaking via 'trans- 
formation acoustics'. The transformation media proposed in |10|l24j are rather singu- 
lar. This poses much challenge to both theoretical analysis and practical fabrication. 
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In order to avoid the singular structures, several regularized approximate cloaking 
schemes are proposed in [6l [12l [131 [13 125] • The basic idea is to introduce regular- 
ization into the singular transformation underlying the ideal cloaking, and instead of 
the perfect invisibility, one would consider the 'near- invisibility' depending on a reg- 
ularization parameter. Our study is closely related to the one introduced in [13] for 
approximate cloaking in EIT, where the 'blow-up-a-point' transformation in \10\ [23] 
is regularized to be the 'blow-up-a-small-region' transformation. The idea was further 
explored in |12| and [T7] for the Helmholtz equation. In [T7], the author imposed a 
homogeneous Dirichlet boundary condition at the inner edge of the cloak and showed 
that the 'blow-up-a-small-region' construction gives successful near-cloak. In |12J, the 
authors introduced a special lossy-layer between the cloaked region and the cloaking 
region, and also showed that the 'blow-up-a-small-region' construction gives success- 
ful near-cloak. For both cloaking constructions, it was shown that the near-cloaks 
come, respectively, within 1/| lnp| in 2D and p in 3D of the perfect cloaking, where p 
is the relative size of the small region being blown-up for the construction and plays 
the role of a regularization parameter. These estimates are also shown to be optimal 
for their constructions. 

It is worth noting that if one lets the loss parameter in [12] go to infinity, this 
limit corresponds to the imposition of a homogeneous Dirichlet boundary condition 
at the inner edge of the cloak. On the other hand, the imposition of a homogeneous 
Dirichlet boundary condition at the inner edge of the cloak is equivalent to employing 
a sound-soft layer right outside the cloaked region. In this sense, the lossy layer lining 
in [l2j is a finite realization of the sound-soft lining in [17]. We would like to emphasize 
that employing some special lining is necessary for the near-cloak construction, since 
otherwise it is shown in [12] that there exists certain resonant inclusion which defies 
any attempt to achieve near-cloak. 

In this work, we shall impose a homogeneous Neumann boundary condition on the 
inner edge of the cloak, which amounts to employing a sound-hard layer right outside 
the cloaked region. The cloaking scheme is referred to as an SH construction. For 
the SH construction, we show that one could achieve significantly enhanced cloaking 
performance. Actually, it is shown that one could achieve, respectively, in 2D 
and p^ in 3D within the perfect cloaking for such construction. We then develop 
a cloaking scheme by making use of a well-chosen lossy layer lining. The properly 
designed lossy-layer could be taken as a finite realization of the sound-hard lining. 
The cloaking scheme is referred to as an FSH construction. The FSH construction is 
shown to possess the same cloaking performance as the SH construction. 

The analysis of cloaking must specify the type of exterior measurements. In 
|6l [T2l I13j . the near-cloaks are assessed in terms of boundary measurement encoded 
into the boundary Dirichlet-to-Neumann (DtN) map. The scattering measurement is 
considered for the near-cloak in [T7]. In the present work, we shall assess our near- 
cloak construction with respect to scattering measurement encoded into the scattering 
amplitude. Nonetheless, by |2H I22| . it can be shown, at least heuristically, that 
knowing the scattering amplitude amounts to knowing the boundary DtN map. 

In this paper, we focus entirely on transformation-optics-approach in constructing 
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cloaking devices. But we would like to mention in passing the other promising cloaking 
schemes including the one based on anomalous localized resonance ^Oj , and another 
one based on special (object-dependent) coatings [Ij. It is also interesting to note 
a recent work in [3], where the authors implement multi-coatings to enhance the 
near-cloak in EXT. 

The rest of the paper is organized as follows. In Section 2, we develop the cloaking 
scheme by employing the sound-hard lining. In Section 3, we present the cloaking 
scheme with a properly designed lossy layer. Section 4 is devoted to discussions on 
different cloaking schemes. In Section 5, we present the numerical examples. 



2 Transformation acoustics and cloaking construction 

Let q G L°^(M^) be a scalar function and a = (o"*-')^-^]^ G Sym(A^) be a symmetric- 
matrix- valued function on M^, which is bounded in the sense that, for some constants 
< Co < Co < oo, 

coc'^e < < coe^e (2.1) 

for all X G and ^ G M^. In acoustics, a^"^ and respectively, represent the 
mass density tensor and the bulk modulus of a regular acoustic medium. We shall 
denote {M''^;ct, (7} an acoustic medium as described above. It is assumed that the 
inhomogeneity of the acoustic medium is compactly supported, namely, a = I and 
g = 1 in M''^\Q with Vl a bounded Lipschitz domain in M^. In M^, the time-harmonic 
acoustic wave propagation is governed by the heterogeneous Helmholtz equation 

div(o-VM) + k'^qu = 0, (2.2) 

where A; > represents the wave number. Stationary scattering theory is to seek a 



solution to (2.2) admitting the following asymptotic development 



= e"-^ + , uTLo 1^ fo> 4 V fo ) ^ ^ (2-3) 



x\{N-i)i2 y y\xy \\X 

where ^ = kd with d G S^^^. A{x,d) with x := x/\x\ is the so-called scattering am- 
plitude. An important problem arising in practical application is to recover {$7; a, q} 
from the measurement of the corresponding scattering amplitude. In the following, 
for clarity and also the convenience of our subsequent study, we give a more de- 
tailed description of the scattering problem. We shall let u^{x) := e^^^''^ denote a 
time-harmonic plane wave, where d G S^~^ denotes the incident direction. Let u*""* 
and u^^^ denote the total wave fields inside and outside the inhomogeneous medium. 
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respectively, which satisfy the following PDE system 



V • (dVu*"*) + k^qu'"* = in n, 



2„,ext 



in 



u \dn 



u Ian, 



,int 



u^^^{x) = u'{x) + u''{x), xe R^\n 

hm r(^-i)/2|^_iA;n4=0, 
or 



diL 



ext 



an 



an 



(2.4) 



where r = \x\ for x G M^. We know e ///^^(M^) (see, e.g. [II1[I9]) and clearly, 
A{x,d) can be read off from the large |x| asymptotics of u'^. 

In this paper, we shall be concerned with the construction of a layer of cloaking 
medium which makes the inside scatterer invisible to scattering amplitude. To that 
end, we present a quick discussion on transformation acoustics. Let x = F(x) : Q — )• 
Q be a bi-Lipschitz and orientation-preserving mapping. For an acoustic medium 
{SI; <T, q}, we let the push-forwarded medium be defined as 



where 



{Q; cr,q} = F^{Q; a, q} := {Q; F^a, F^q}, 



a{x) = F,a{x) := jMa{x)M^\^=p-^^) 
q{x) = F^q{x) := q{x)/ J\^^p-i^^^ 



(2.5) 



(2.6) 



and M = {dxi/dxj)^j^i, J = det(M). Then u G H^{Q) solves the Helmholtz equa- 
tion 

V • (crVu) + k'^qu = on r2, 
if and only if the pull-back field u = {F~^)*u := uo F^^ G H^{Q) solves 

V • (aVu) + k'^qu = 0. 

We have made use of V and V to distinguish the differentiations respectively in x- 
and x-coordinates. We refer to [12^ [T7] for a proof of this invariance. 

We are in a position to construct the cloaking device. In the sequel, \et D <£ 0, 
be a Lipschitz domain such that 0,\D is connected. W.L.O.G., we assume that D 
contains the origin. Let p > be sufficiently small and Dp := {px; x G D}. Suppose 



Fp : n\Dp n\D, 



(2.7) 



which is a bi-Lipschitz and orientation-preserving mapping, and Fplg^ = Identity. A 
celebrated example of such blow-up mapping is given by 



y = Fp{x) 



R2 + — \x\ 



R2-P 



R2-P 



X 

\x\ 



p < Ri < R2 



(2. 
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which blows-up the central bah Bp to Buj^ within B^.^ ■ Now, we set 



(2.9) 



Let M D represent the cloaked region. Then we claim the following construction 
gives a near-cloaking device 



a sound-hard layer 
arbitrary target object 



in M^\0; 
in n\D- 
in D\M; 
in M. 



(2.10) 



In (2.10), by a sound-hard layer we mean a layer of material which prevents acoustic 



wave from penetrating inside and the normal velocity of the underlying wave field 
vanishes on the exterior boundary of the layer. The wave equation governing the 



wave scattering corresponding to the cloaking device constructed in (2.10) is 



in 



V • (crVn) + k'^qu 

N 

^{aPy^VidjU = Q on dD, 



(2.11) 



where v = is the exterior unit normal vector to dD. In the following, we shall 

let A{x,d) denote the scattering amplitude to the PDE system (2.11) corresponding 
to the cloaking device. Let 

F = Fp on 9\Dp- Identity on R^\9.. 

Set V = F*u G Hl^^(R^\Dp) and v^{x) := v{x) — e'^^''^. By transformation acoustics, 
together with straightforward calculations, one can show that 



(A + P)v = in 
=0, 

v'^{x) + e*^'^ X G 
— ikv 

r^oo I or 



dv 

lim r(^-iV2 



0. 



(2.12) 



In terms of the terminologies in [T7], (2.11) describes the scattering in the physical 
space and (2.12) describes the scattering in the virtual space. Since v = urn. M^\0, 



we see the scattering in the physical space is the same as that in the virtual space. 
That is, A{x, d) could also be read off from the large |x| asymptotics of v*. Next, we 
give one of the main results of this paper, which justifies the near-invisibility of the 
above construction. 
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Theorem 2.1. There exists po > such that when p < pQ the scattering amplitude 
A{x, d) to (2.1^ satisfies 



\A{x,d)\<Cp^, x,d£S 



(2.13) 



where C is a constant dependent only on k, pQ and D, hut completely independent of 
P- 

Proof. In order to ease the exposition, we shall only prove the theorem for = 2, 3. 
But we would like to emphasize that for > 3, the proof follows by completely 
similar arguments. 
Let 

. / , \{N-2)/2 
W K \ ^(1) 



G{x) 



4 V 27r x 



(Af-2)/2 



{k\x\) 



(2.14) 



be the outgoing Green's function. By Green's representation, we know for x G M^\D 



v'{x) 



\ dG{x-y) ^.. dv'jy) 
--{ICvn{x)+g{x), 



ds{y) 



(2.15) 



where we have set 



{ICv^ix) :- 
g{x) :-- 



dDo 



dv{y) 
dD, dv{y) 



(2.16) 
(2.17) 



Clearly, v%x)\aD, e H^^dDp). By the jump properties of the double-layer potential 
operator /C (cf. [19]), we have from (2.15) 



'-v'{x) = iJCv'){x) + g{x), X £ dDp 



(2.18) 



Let x' = x/p, then (2.18) is read as 



y{px') = {lCv'){px') + g{px'), x' G dD. 



(2.19) 
(2.20) 



Next, we claim 

hip ■)\\L^{dD) < Cp, 

where C remains uniform as p — >• 0+. In the sequel, we shall make use of the following 
asymptotic developments of the 2D G{x) (cf. [5]), 



1 i 

Gix) = In Ixl + - 

^ ' 27r ' ' 4 27r 2 27r 



1 X^^_^ + 0{\x\Hn\x\) 



(2.21) 
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for — )• 0, where E is the Euler's constant. Moreover, we know that 



ik\x\ 

G(x) = when N = 3. 



(2.22) 



In order to prove (2.20), we first assume x G dDtp with 1 < t < 2. Then by Green's 
formula, we have 

Q^iky-d 

G{x - y) r. . . ds{y) 

dD, My) 

: [ Aye'''y'^G{x -y)dy+ [ VyG{x - y) ■ Vye'^y^ dy (2.23) 
Jd„ Jd„ 



e'ky'Gix -y)dy+ I VyG{x - y) ■ Vye'^y" dy. 



Da 



iky-d 



Do 



Using (2.21) and (2.22), we have for x' £ Dt 



\g,{px')\=k^\ / e'''y''G{px'-y)dy\ 
<k^ [ \Gipix'-y'))p^dy'\. 



D 



By the mapping properties of volume potential operator, one has 

\9iix)\c{dDtp) < Cp, 
where C is independent of p and t. In like manner, one can show that 

\92ix)\cidD,,) = I / VyG{x - y) ■ Vye'^y" dy\c^aD,,) < Cp. 
Jdo 



By pj24^ and i\L2E\ i, we see 

\9{x)\c{Dtp) < cp. 

Next, by the mapping property of single-layer potential operator, we know 

9{x)\9D, = ^hm {g{x)\QD,p) , 



which together with (2.26) implies (2.20). 



(2.24) 

(2.25) 
(2.26) 
(2.27) 



We proceed to the integral equation (2.18). First, by using change of variables in 
integrals, it is straightforward to show that 

ilCv'){px') = ilCov'ip-)){px') + {Rv%p-)){px'), x'GdD, (2.28) 

where /Cq is an integral operator with the kernel given by 



Go{x-y) = < 



In |x — w| N = 2, 

2. I 

^ 47r |a; — yl 
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which is the fundamental solution to —A; and R satisfies 

\\R\\c{L'2{D),L2{D)) 



< 



plnp when N = 2; 
p when = 3. 



(2.29) 



Hence, the integral equation (|2.19|) can be reformulated as 

{px') = g{px) X G D. 



{h-)Co-R)v'{p 



(2.30) 



By the well-known result in [27], I — |/Co is invertible from L^{dD) to LF'{dD). Then, 
by using (2.20), we have from (2.30) that 

\\v'{p ■)\\mdD) < CMp ■)\\LHaD) < Cp. (2.31) 

Noting ||t'^(-)llL2(Q£)p) = /3(^"^^/^||?j*(p •)llL2(a£)), we further have from (2.31) that 



\V \\L^(dDp) 



< 



r C7p3/2 ^ ^ 2, 
ycp^ TV = 3. 



(2.32) 



Finally, by letting \x\ — >• oo in (2.15), we have 
A{x, d) = 7 



dDp 



fip—ikx-y fipiky-d 



dv{y) 



dsiy), (2.33) 



where 7 = e'4 /y/Sirk when N = 2, and 7 = l/An when = 3. By using (2.32) and 



the Schwarz inequality in (|2.33), we have 



dDp 



dv{y) 



v%y) ds{y) 



<Cp^. 



(2.34) 



By Green's formula, we have 



-ikx-y 



dDp 



Q^iky-d 

dv{y) 



ds{y) 



[ (^Ae'''y'^)e-'''^-y dy+ [ Ve'^y^ ■ Ve'^'^^-y dy 

JDp JDp 



(2.35) 



=\k^{d-x-l) / e'''^'^~^^-y dy\ 

JDp 

<Cp''. 



By ( |2.33D , \2M\ and ( |2.35D , we have \2.l^ . 
The proof is completed. 



□ 
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By Theorem 2.1, we know the construction (2.10) gives a near-invisibihty cloaking 
within of the perfect cloaking. For a special case by taking Dp = Bp, namely, the 
central ball of radius /c > 0, using wave functions expansion, one can show (cf. |18j ) 



A{x,d) = -e"*4 . 

in M^, where 9 = d); and 
A{x,d) -- 




2 



+ 2f ^^^cosn^ 



n=0 



h^n^ {kp) 



cos ( 



(2.36) 



(2.37) 



in M^, where P„ is the Legendre polynomial of degree n. By the asymptotic develop- 



ments of spherical Bessel functions (cf. [18j), one has from (2.37) 

{kpf + 0{{kpf). 



A{x,d) 



k 



cos ( 



Similarly, by (2.36) one can show that in 



A{x, d) 



2lT ( cos( 

k 



{kpf + 0{{kpf). 



(2.38) 



(2.39) 



By (2.38) and (2.39), it is readily seen that the estimates in Theorem 2.1 are optimal 
for full scattering measurements, namely, x G S^~^ and d G Nevertheless, it 

is interesting to note from (2.38) and ( |2.39 ) that for some specific scattering mea- 
surements, e.g., A{x,d) with Z{x,d) = ibarccos| in 3D, and with Z{x,d) = ib| 
in 2D, one would have even more enhanced invisibility cloaking effects. Physically 
speaking, the cloaking effect would be stronger in the backward scattering region, i.e. 
|Z(x,(i)| > tt/2, than that in the forward scattering region, i.e. |Z(x,d)| > tt/2. This 



could be partly seen from (2.33)-(2.35), and (2.38)-(2.39). 



3 Near-cloak construction with a lossy layer 

In this section, we shall develop a lossy approximate cloaking scheme. To that end, 
we first introduce the following transformation T : M.^ — ?■ M^, 



y = T{x) :-- 



where Fp is given in (2.7). Let 



{M^;a,g} 



X 

p 



for X G R^\n, 
for X £ ^\Dp, 
for X G Dp, 



LI 



in M.'^\n, 
in n\D, 
T^ai,T^qi in D\Di/2, 



in Dii2, 



(3.1) 



(3.2) 
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be the cloaking device in the physical space. Here, {Di/2',cr'^,qa} is an arbitrary but 
regular medium, which represents the target object being cloaked; and 

{Dp\Dp/2;^i,qi} (3.3) 

is a lossy layer whose parameters shall be specified in the following. Similar to our 
earlier argument in Section 3, by transformation acoustics we see that the scattering 
amplitude in the physical space corresponding to the cloaking device is the same as 
the one in the virtual space. In the virtual space, the wave scattering is governed by 
the following PDE system 



■(A + A;2)u = in R^\D^ 

V ■{aiVu)+k'^qiu = in Dp\Dp/2, (3.4) 

V • (fJaVu) + k'^QaU = in D 



V2' 



{Dp/2;aa,qa} = {T-')ADy2; q'a) (3-5) 



where 

is arbitrary but regular; and u G i?;^^^(M^) satisfies 



We shall choose 



u{x) = e^'^^-'^ + u'{x) for x G M^\D 
lim A^-^)/^!^-iku'[=0. 

r— >oo I or 



ai = Cp'^+'^^I and qi = a + ib (3.6) 



with C, 5, a, b some fixed positive constants in our construction. We shall show that 
it will yield an enhanced approximate cloaking scheme. However, the proof for the 
general case with general geometry and arbitrary cloaked contents is rather technical 
and lengthy, which we choose to extend to full details in our forthcoming paper |15| . 
In this section, we shall consider the special case with spherical geometry and uniform 
cloaked contents. In the sequel, we let Dp = Bp and, a'^ and q'^ be arbitrary positive 
constants but independent of p. With a bit abusing of notation, we shall also write 



ai = Cp 



2+2S 



By (3.5), one can show by direct calculations that in the virtual space 



O'a, qa 



in Dp/2 when N = 2 



p- 

pN 



in Dp/2 when 



For the PDE system (3.4), we let u = uq in \Dp, u = U2 in Dp\Dp/2 and 



Uint in Dp/2. By transmission conditions, we have 



du2 duo on /Qv^ 
U2 = uo and ai^— = —- on dDp, (3.7) 

OU Of 
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and 



du2 dUint 
= Uint and ai— — = aa^^^ — on oD 



dv dv 



p/2- 



(3.8) 



Set k = kx — and h 



.:2 = K^i^. We choose the complex branch of k such that 



%k) > 0. 

We first consider the 2D case. By [3], one has the following series expansions 



n=—oo n=—QO 
oo oo 

U2{x)= J2 a„J„(^|x|)e*"^+ J2 bnHW{k\x\)e 

n=—oo 
oo 

Uint = Y CnJn{k2\x\)e 



ind 



And 



(3.9) 



„in6 



By (3.7) and (3.9), we have 

janJn{kp) + hnH^\kp) = i^Jnikp) + dnH^n\kp) 



ai[ankJ'^{kp) + bnkH^n' {kp)] = t^kJ'^ikp) + dnkH^^^' (kp). 



(3.10) 



In like manner, by (3.8) and (3.9) we have 



anJn{kp/2) + hnHn{l){kp/2) = CnJn{k2P/2) 

ai[ankJn{kp/2) + hnkH^^^' Ckp/2)] = aaCnk2Jn{k2p/2). 



(3.11) 



Here J'^{kp) = ^1,=^, and j;(fcp/2) = ^L^^,/^. H)^' {kp), H\^' {kp/2), 
J'^{kp), Hn (kp), J'^{k2p/2), Hn {k2p/2) are understood in the s ame s ense. By 



r(i)' 



letting Co = l/^/alqi and A = ^qa^a = k2 = ^, (|3.10|) and (|3.1l|) are 



read as 



and 



anJn{kp) + bnHi^\kp) = i^Jnikp) + dnHl^\kp) 



anJ'nikp) + bnHi'^ (kp)] = Co[i" ^(A;^) + dnH^'' (kp)] 

(anJn{kp/2) + hnH^n\kp/2) = CnJn{k2p/2) 
\[anJ'n{kp/2) + bnH^^^' Ckp/2)] = CoACnJ'n{k2p/2). 



(3.12) 



(3.13) 



Noting k2p = k^/q^Ja^, if Jn{k2p/2) ^ then by direct calculations, we have from 
( 3l3l > 

anJnCkp/2) + bnHi^\kp/2) 



and 



Jnik2p/2) 

J'^{kp/2)-CoA'-f;^^Ukp/2) 
' R'^Hkpl2)-C,A'j^^H^^\kpl2) 



a„,. 



(3.14) 
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In case Jn{k2p/2) = 0, we have from the first equation in (3.13) that 

J„(^p/2) 



bn 



Hi'\~kp/2) 



(3.15) 



Let To denote the fraction in (3.14) or (3.15) and hence bn = Tgan- Plugging bn into 
(3.12), we have by straightforward calculations 



and 



dn 



_ i'^Jn{kp)+dnH'n '{kp) 

Jnikp) + ToHL^\kp) 

i^J'n{kp) - 1 j;(M+Tfogi^''(M jnj (fe ) 

Hi'^ikp) - 1 ^4(M+ToHi^)'(M ff(i)(fe )• 



(3.16) 



We next investigate the asymptotic development of 



n{(Ti,p) :-- 



1 j;(M + ToM')'(M 

C^o Jn{kp) + ToH^^\kp) 



1 j;(fcp) 1 gi'^'(fcp) 

Co J„(fcp) Co ° J„(fcp) 
^ ^° J„(fcp) 



(3.17) 



/ (1) 

Clearly, we only need study the asymptotic behaviors of j- '\^%^\ , -^Tq^^ytt^ 
'^o ^T ■ We note the following fact due to (3.6), 

^{kp) — )• +00 and Q{kp) — )■ +oo as p — )• 0"*", 

which implies the following asymptotic developments for z = kp (see formulas 9.2.1, 
9.2.3 in [2]), 



nvr 
2 



f) + el^(-)IO(|z|-i), 



arg z\ < TT 



riTT TT \ 

2 i', — TT < argz < 27r 



(3.18) 



We also note the following recurrence relation (see formula 9.1.27 in [2]) for the 
subsequent use, 

(3.19) 



n 



B'Jz) = -Bn{z) - Bn+l{z) 



where ;S„(z) = Jn{z) or Hn\z). Since Jn{z 



-l)"J„(z) andF, 



(see formula 9.1.5 in [2j), we onl y need consider the case with n > in the sequel. 
Noting cos(z) = — , by (3.18) we further have as ^{z),'is{z) — c« 



Jniz) ~ 
Hi^'\z) ' 

\hI^^'{z) 



2 



arg z\ < IT 
-TT < arg z < 2n 
TT < arg z < 2tt 



(3.20) 
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It is rem arked that the third equation in (3.20) is obtained by using the recurrence 
relation (3.19). Hence we see for z = fcp, Jn{,z) blows up exponentially, while Hn\z) 
decreases exponentially as /o — )• 0+. In the sequel, we let \/a + ib = a + i/3 with 



a,j3>0. For 



1 J4(fcp) 
Co J„(fep) ■ 



as /9 — )• 0+, by (3.20) we have 



1 _ 1 f/nCkp) - Jn+l{kp) 

Co Jn(kp) C'O Jn{kp) 



Co 



For To ^^^7^^, if To is the fraction in (|3.14|), then we have 



(3.21) 



T, 



H^r^\kp) ^ Ukpll) - CoA'-f^^Ukp/2) Hi'\~kp) 
' JnCkp) ^«'(fcp/2)-Co^J(ggi?W(^p/2) JnCkp) 

J'^{kp/2) - CoAj^^Ukp/2) Hi'\kp) 



Jn{kp) 



H^rthhm-coAif^^-^^^ 



■.=yi X 3^2. 



m:'{kp/2) 

(3.22) 



By straightforward asymptotic analysis, one can show as p — ?• 0"* 

1 4(fc2P/2) 



13^1 



j;(W2) - Co^^igf^ J.(W2) 



Jn{kp) 



>/2 



CoA 



J: 



l{k2p/2) \ Ukp/2) 

^Jn{k2p/2)) J^{kp) 
<C(1 + (71 + 1)10 + ^/31/5-2-5 



(3.23) 



where C is a constant independent of p. That is, |3^i| decreases to more quickly 
than p"^ for any r > 0. Similarly, one can show that |3^2| decreases to more quickly 



than p^ for any r > 0. Furthermore, if Tq is the fraction in (3.15), by completely 
similar arguments one can show that Tq " ^ ^' also decays more quickly than for 
any r > 0. Hence, by ( 3.17 )-( 3.23) and ct; = p^^"^^ , we have 



^(o-/,p) 



1 j;(M 



1 



e*^/2(c, + i/3)pi+<5 as /3^0+. (3.24) 



Co Jn{kp) Co 
Now, by ( |3.24D and ( |3.16[ ) we have 

i'^J'nikp) + e'^/2(a + /3i)p^+'^i" Jn(M 
M^^'(fcp) + e^-/2(a + Pi)p^+^H^n\kp) 



dn 



(3.25) 
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Using the asymptotic developments of Bessel functions and their derivatives (cf. jl8j). 
we further have 



1+5 



do ~ 



-kp/2 + e''^/^{a + il3)p 



n = 0, 



irkp 
i"n{kp) 
2"r 



(n+1) ''Plt^ * 2"r{n+l) 



(3.26) 



, n G N, 



which imply 



and 



do - 




dn - 







(io~0(/o^), n = 
dn^-Oip^"^), nGN 



< 7rA;|a + f/3|/92+'^, 
47r|a + i/3| 1 



(3.27) 



(3.28) 



< 



Since 



^1+5 (2nn!)2 

oo 

u'{x)= dnH^^\k\x\)e 



{kp) 



2n+2+S 



ind 



(3.29) 



by taking \x\ — )• +oo, together with (|3.27|), one has by direct calculations that the 

(3.30) 



corresponding scattering amplitude satisfies 

\A{x,d)\ < Cp^, 



where C is a positive constant that remains uniform as p — > O"*". That is, the con- 
struction (3.2) gives an near-cloaking device within p^ of the ideal cloaking. 

Now, we look into the series representation of the scattered wave filed (3.29). If 
Dp is a sound-hard obstacle, the scattered wave corresponding to e^^^''^ is given by 
(see eqn. (3.19) in [IS]) 



tt*(x) 



Hi'^'ikp) 



H^^\kx)e 



ind 



(3.31) 



By (3.28 ), (3.29) and (3.31 ), one has by direct verifications that for p sufficiently small 

\u%x)-u'{x)\<Cp^+^ (3.32) 

for any x G ^\B^f^ with > p a fixed constant, where C depends only on k but 
independent of p. We remark that the estimate in (3.32) would reduce to Cp^^^ 
for X £ M.'^\Dp. Hence, the construction (3.2) is actually a finite realization of the 
sound-hard lining construction (2.10). 
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The 3D case could be proved following similar arguments, which we shall sketch 
in the following. We shall make use of the same notations uq, U2, Uint, k and /c2 etc.. 
We have the following series representations of the wave fields, 

oo n CO m 

=E E i''^^Wid)jn{k\x\)Y;ri^) + Y^ d^hi'Hk\x\)Y;ri£), 



n=0 m=—n 
oo n 



n=0 n=—m 



oo n 



n=Q m=—n n=0 m=—n 

oo n 

u^nt{x)=Y, E CCik2\x\)YriS:). 

n=0 m=—n 

By the transmission conditions, respectively on dDp and dDp/2y we have 
<jnCkp) + Kh^n\kp) = i^'^TTYTid^jnikp) + d^h^^\kp), 



(3.33) 



and 



kal'j'^{kp) + kh'^h^n^' {kp)\ = i''AnkY^{d)3'n{kp) + kd^^h^n^' {kp), 
'a^Ukp/2) + h^h\i\kp/2) = CUk2p/2), 



ka^j'^(kp/2) + kh^h^^^ {kp/2) \ =aaCk2Uk2p/2). 



(3.34) 



(3.35) 



For this 3D case, we have A = y^cr^g^/p^, ^2 = ^ y fr- Similar to the 2D case (see 
(3.13), (3.14) and (3.15)), we would solve the linear systems ( |3.34 ) and (3.35) and we 
need to distinguish two cases jn{k2p/2) / and jn{k2p/2) = 0. In the following, we 
only present the more complicated case with jn{k2p/2) ^ and the other case could 
be handled in a completely similar manner. By solving (3.34) and (3.35), we have 



hW\kp) - 1 .;(M+To.i;> (M ;,(l)(fc ) 

f^Ckp/2) - CoAf^^Ukpm 
h^^y(kp/2)-CoA&^,h'i\kp/2) 



where 



Tn : 



jn{k2p/2) 

By similar asymptotic analyses as those for the 2D case, one can show that as p — )• O"*", 



1 j'^Ckp) + '^oh'n^'Ckp) 

Co jn{kp) + '^0hl^\kp) 



1 
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Then by ( 3.36[ ) we have 



4- 



j(2n)!(n+l) „i7T/2(^. i AO) /Zf: ^IM^ 



2"n!(A:p)" 
{2n+l)! 



(3.37) 



which in turn imphes 



and 



(3.38) 



d'. 



47rY,^{d)j',{kp) 



< 8TTk'^\a + i(3\\Y^{d)\p^+^, 



h^n^ {hp) 



< 



4(2"n!)247r|y7((i)||a + i/3|/c- 
(2n)!(2n + l)!(n + l) 



-1-5 



-(M 



2n+3+5 



(3.39) 



With the above preparations, one can show that the corresponding scattering amph- 
tude corresponding to the cloaking device (3.2) in satisfies 



\A{x,d)\ < Cp^ 



(3.40) 



where C remains uniform as /c — )• 0^. That is, (3.2) gives a near-cloaking device 
within p^ of the ideal cloaking. Furthermore, by comparing the scattered wave field 
in (3.33) to that of a 3D sound-hard ball Bp 



[18]), together with (3.39), one 
can show that the deviation between them is within p^~^^ in M.^\Beo for any fixed 



eo > p, and within \Bp. Hence again, we come to the conclusion that the 



construction (3.2) is a finite realization of the sound-hard lining construction (2.10). 



In summary, we have shown in this section that 



- Bp and {Dp\Dp/2;ai,qi} be given by (3.6), and {Dp/2; 



Theorem 3.1. Let 

Qa} be arbitrary but uniform. Then the construction (3.2) produces a near-cloaking 
device within p^ of ideal cloaking. Furthermore, (3.2) is a finite realization of (2.10) 



in the sense that the scattered wave fields corresponding to (3.2) and (2.10), respec- 
tively, deviate within p^~^^ outside the cloaking device. 



As we remarked earlier that Theorem 3.1 equally holds for the general case with 



general geometry and arbitrary cloaked contents (see |15j). Actually, in the subse- 
quent section, our numerical examples are presented for variable cloaked contents. 
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4 Some discussion on different cloaking schemes 



As we discussed earlier in Section [T| two different cloaking schemes were developed in 
|12j and |17] . One of the key ingredients is to introduce a special layer between the 
cloaked region and the cloaking region, which are respectively -D1/2 and i}\D in the 



physical space described in (2.10). They correspond, respectively, -Dp/2 and ^l\Dp in 
the virtual space. In [12], the authors introduce a lossy layer occupying D\Di/2 and 
they showed that the lossy layer is indispensable to achieve successful near-cloak. In 
virtual space, the lossy layer in [T2] is given as 

{Dp\Dp/2;I,l + if3}, P^p-^. (4.1) 

If one lets the loss parameter (3 go to infinity, the limit corresponds to the imposition 
of a homogeneous Dirichlet boundary condition, which is the one considered in [T7] . It 
is interesting to mention that the improvement of cloaking performance by imposing 
specific boundary condition on the cloaking interface is also considered in [9]. In this 
context, imposition of a homogeneous Dirichlet boundary condition amounts to the 



lining of a layer of sound-soft material in Dp\Dpi2- In this sense, the lossy layer (4.1 ) 
is a finite realization of the sound-soft layer lining. We shall refer the construction 
in ^1 with a sound-soft layer as SS construction, and the one in [T7] with a finite 
realization of sound-soft layer as FSS construction. 



In our cloaking construction (2.10), the inclusion of a sound-hard layer will sig- 
nificantly enhance the cloaking performance, and as specified in the Introduction we 
call this scheme an SH construction. In order to achieve a finite realization of the 
sound-hard layer, we utilize the following lossy layer 

{Dp\Dp/2;(Jua + ih] , ai^p^+^I, a ~ 1,6-1. (4.2) 



Clearly, the lossy layer (4.2) is of significant physical and mathematical nature from 



(4.1), and it could produce significantly improved near-cloaking performances. This 



losses scheme is called an FSH construction. 



5 Numerical examples 

In this section, we present some numerical experiments to demonstrate the theoret- 
ical results established in the previous sections. There are totally four near-cloak 
schemes investigated, namely SS (sound-soft lining), FSS (finite sound-soft lining), 
SH (sound-hard lining) and FSH (finite sound- hard lining). 

First, all the numerical experiments are conducted in and we choose D and Dp 
as -Bri and Bp, respectively. Some key parameters are set as follows: i?i = 2, i?2 = 3 
and i?3 = 4, thus Rq = i?i/2 = 1. Experimental settings are shown in Figure [Tj 
For the setting with respect to the scattering measurement, a Cartesian PML layer 
of width 1 is attached to the square (—5,5)^ to truncate the whole space into a 
finite domain with scattering boundary conditions enforced on the outer boundary. 
(—5,5)^ \ -B/jg and Bji^ \ Bji^ are homogeneous surrounding media. Br^ \ Br^ is 
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Figure 1: Experimental settings for scattering measurement tests. Left: Tests with 
SS and SH lining. Right: Tests with FSS and FSH lining. 



the cloaking medium. Br^ \ Br^, Bpi^ are the lossy layer and the cloked region, 
respectively. The regularization parameter p ranges from 0.5 to 10~^. We fix the 
wave number k = 2, and the incident direction d = (1,0)^. 

In addition, for Scheme FSS, the lossy parameters are chosen according to |12] . 
i.e., ai = I, qi = + 2.5/j~^i) in Br^ XBr-^. For Scheme FSH, the lossy parameters 
are chosen by (3.2), (3.6) and the definition of the map (2.8), namely ai = p^'^I, 
qi = /9^(3 + 2i) in the cloaking medium of the physical space, namely C = 1, 6 = 0.5, 
a = 3 and 6 = 2. 

For Schemes FSS and FSH, We choose the medium coefficients in the cloaked 
region to be variable functions da = 1 + 0.3 x^y^, qa = 5 + x. For Scheme SH, the 
incident wave is shifted to the right by a tenth of the wave length such that the origin 
is not in the valley nor peak and the incident wave does not vanishes at the origin. 



From the transformation map F defined in (2.8), we can determine the cloaking 



medium coefficients by (3.2) for any p < Ri. It is emphasized that as p tends to 
zero, the singularity of a increases very fast. For example, it is observed that an 
grows asymptotically as 1/p and exhibits a thin layer of large values except some 
constant background as p decreases, which requires many local refinements for better 
resolution and accounts for the boundary layer in the finite element solutions. 

The scattering amplitude, or the far field data, can be computed in the following 
way (cf. [5]). We can represent the scattered wave in the boundary integral form 

d^{x, y) ^, du'iy) 2 



u'{x)= u%y) -^x,y)—^ds{y), xeR'\BR^. (5.1) 



9Bfl3 du{y) ' (9i/(y) 



Making use of the asymptotic expansion of the fundamental solution $(x,?/) 



iHQ^\k\x — y|)/4, we then have 



exp(-ivr/4) f ^,,,^,,9e ^^^-y ^_,kx-y du' (y) 



A{x) := " ^ ' / : , , - e-*'=^-2/r^l^ ds{y), (5.2) 

^ ^ VskTT JdBn, My) My) 

where x = x/\x\ is the observation direction and ^{y) is the outward unit normal 
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vector pointing to the infinity. That is, the far field data can be approximated by the 
numerical quadrature of (5.2) using the Cauchy data of the scattered wave on the 
boundary of B3. 

Two groups of tests are carried out in the following. 



5.1 Near-cloaks of scheme SS and scheme FSS 

It is pointed out that as p tends to zero, the medium coefficients exhibit more and 
more singularity close to dBji^ . Thus finite element solutions require more degrees of 
freedom to resolve the boundary layer close to . Either the sound-soft boundary 
condition or the FSS lossy layer enforces a strong influence to the local behavior of 
the transmitted wave close to dBji^ , which amounts to introducing a point source in 
the virtual space. Such virtual point source makes the boundary data show a very 
slow convergence as we decrease p. To tackle with the difficulty of the boundary layer 
amounts to more than two millions of DOF's to achieve the desired relative error 
tolerance when p = 10~^. The same observation holds for the other schemes. To 
avoid the influence of the finite element discretization error, all the following tests are 
carried out on the finest mesh to show the sharpness of the theoretical error bounds 
in terms of p. 

The scattered wave and the transmitted wave are plotted for Schemes SS and 
FSS, respectively, for different p's in Figures [2j 




^1 






Figure 2: The scattered wave and the transmitted wave (real part) with respect 
to p = 10~^ and 10~^ from top to bottom, respectively. Left: Scheme SS; Right: 
Scheme FSS. 



The scattering measurement data are generated according to (5.2) by solving the 



Helmholtz system on the truncated domain with a PML layer. In the discrete sense, 
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i/i;osio(p)i 



l/l'oSioWl 



Figure 3: Convergence history of scattering measurement data versus 1/| log^^Q p\ for 
Scheme SS (Left) and Scheme FSS (Right). 



all the norms are equivalent. So the discrete maximum norm are used to measure 
the decay rate of the far field data A{x), namely the maximum of the modulus of 
A{x) are taken over 100 equidistant observation direction on S^. We investigate the 
convergence rate by testing p = 1/2-^, j = 1,2, ... ,7 for Schemes SS and FSS. By 
plotting in Figure [I the convergence history of the discrete maximum norm of A{x) 
over with respect to the upper bound 1/| log^oPl) see clearly that Schemes SS 
and FSS schemes indeed achieve near-cloak, but we see that the convergence curve 
demonstrates linearity asymptotically except the first few outliers to the left of the 
plots in Figure [3| which verifies the sharpness of the upper bound established in [12j. 
But the convergence slows down significantly as p decreases, which is reflected by 
the clustering of the blue star points in the curve. More importantly, it can be seen 
from Figure [3] that the difference between the discrete maximum norms for Schemes 
SS and FSS gets smaller as p decreases because the FSS lossy layer tends effectively 
to the sound-soft boundary condition for small p and thus Scheme FSS approaches 
Scheme SS in the limit sense as p — )• 0. 

5.2 Near-cloaks of scheme SH and scheme FSH 

Contrary to Schemes SS and FSS, our new schemes by SH and FSH lining can achieve 
significantly better near-cloak performance. 

Figure |4] shows the scattered wave and the transmitted wave for Scheme SH for 
different p's. Compared with Scheme SS, the scattered wave of Scheme SH decays 
significantly as p decreases. Moreover, the transmitted wave approaches more and 
more like a deformed plane incident wave, the interior value of the transmitted wave 
tends to a constant as /9 — )• 0. 

We plot the scattered wave and the transmitted wave for Scheme FSH for different 
p's in Figure [5] We can observe a thin red layer near the outer boundary of the lossy 
layer, on which both the acoustic potential and the normal flux matches with those 
from the cloaking medium. Compared with Scheme FSS, the scattered wave of Scheme 
FSH decays significantly as p decreases. The transmitted wave approaches more and 
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Figure 4: The scattered wave (real part) and the transmitted wave n* with respect 
to p = 10~^ and 10^^ from top to bottom, respectively, for Scheme SH. 

more like a deformed plane incident wave, and the interior acoustic potential of the 
transmitted wave within the cloaked region vanishes as p — >• 0. In other words, the 
medium scatterer in the cloaked region is approximately isolated from the exterior 
world and has negligible affect on the scattering measurement, which means that we 
achieve significant nearly-invisibility by SH and FSH constructions. 

Finally, we study the convergence history of the discrete maximum norm of A{x) 
with respect to the regularization parameter p. From Figure [6| we see clearly second 
order decay rate of the discrete maximum norm of ^(x) in terms of p for the near-cloak 
construction Schemes SH and FSH, compared with the red reference line of second 
order decay in Figure [6] which confirms the sharpness of our theoretical upper bounds 
for Schemes SH and FSH. Similar to their sound-soft counterparts, it can be seen from 
Figure [6] that the discrete maximum norms for Schemes SH and FSH have nearly 
the same values as p decreases because the FSH lossy layer tends effectively to the 
sound-hard boundary condition for small p and thus Scheme FSH approaches Scheme 
SH in the limit sense as p — >• 0. The significantly improved cloaking performance for 
Schemes SH and FSH makes it easier for engineers to design practical near-cloak 
devices in a variety of industrial applications. 
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